Mechanism of delayed double ionization in a strong laser field 
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When intense laser pulses release correlated electrons, the time delay between the ionizations 
may last more than one laser cycle. We show that this "Recollision-Excitation with Subsequent 
Ionization" pathway originates from the inner electron being promoted to a sticky region by a 
recollision where it is trapped for a long time before ionizing. We identify the mechanism which 
regulates this region, and predict oscillations in the double ionization yield with laser intensity. 

PACS numbers: 32.80.Rm, 05.45.Ac, 32.80.Fb 
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Atoms in strong laser pulses lose their electrons 
through two ionization channels [2]: A sequential one 
(SDI) and its non-sequential counterpart (NSDI). The 
SDI mechanism consists of successive and independent 
removals of the electrons. The recollision (or three-step) 
scenario [3, 4] in which a pre-ionized electron returns to 
the ion core to dislodge a bound electron is the character- 
istic mechanism for NSDI. It turns out that there is a rich 
variety of pathways among NSDI processes, also. At first 
sight, these pathways can be distinguished by the time 
it takes the second electron to ionize. The common vari- 
ant involves little, if any, delay between the recollision 
and ionization. However, this so-called "direct impact 
ionization" [5] is often accompanied by an alternative 
(and less straightforward) road to NSDI called Recolli- 
sion Excitation with Subsequent Ionization (or RESI for 
short |5l4l3J). Recent experiments have shown that RESI 
can be the dominant channel for very short pulses [14]. 
The mechanism for RESI is often attributed to the rec- 
ollision which excites the parent ion, later ionized by the 
laser field with a delay sometimes lasting longer than one 
laser cycle after the recollision and thereby imitating an 
uncorrelated ionization process. 

Here we report results showing that RESI is a purely 
inner electron process and that the recollision, while help- 
ful for increasing its probability, merely enhances what 
is already there, namely the ionization process due to 
extensive state mixing. The reason why the mechanism 
responsible for RESI remains an open question is due to 
the spectator role often erroneously assigned to the inner 
electron. Far from waiting passively for the pre-ionized 
electron to recollide, the inner electron interacts with the 
strong field also, and a true nonsequential description of 
the process reflects the excitation of both electrons. Our 
view, to be detailed below, is summarized by the cartoon 
in Fig. [TJ Hitting the pool of available initial states, a 
recollision may send the inner electron into the ioniza- 
tion continuum directly, thereby causing two electrons 
to emerge with very little time delay (direct impact ion- 
ization) or the recollision may drive the bound electron 
towards the parts of phase space which funnel electrons 
into the ionization continuum - which, given sufficient 



1 electron 

2 dimensiolis 




I 



2.5 


o 




t^ 








N 




^ 


? 


O 




t3 




^ 




3 


1.5 


CD 




^ 




ir> 




rD 




^ 


1 


n 




^ 



0.5^ 



-2 



X (a.u.) 



direct impact ionization 




2 electrons 
2 dimensions 



1.5 2 2.5 3 

Time (laser cycles" 



5 5.5 6 6.5 

Time (laser cycles) 



FIG. 1: Upper panel: Ionization time for Hamiltonian (|2]) 
with two dimensions for / = 3 x 10^^ W • cm~^, 00 = 7r/2 
and 780 nm. Initial conditions are chosen at the energy of 
the Stark saddle [i[ with Px < and Py = 0. After a rec- 
ollision (full line arrows), the inner electron can be directly 
ionized or promoted to an excited state that is ionized with a 
delay (dotted arrow). Lower panels: Typical two-electron tra- 
jectories for direct impact ionization (left) and RESI (right). 
Continuous (dashed) curves denote x (y) coordinates. 



time, the electrons may have reached without the help of 
a recollision. 

This common-sense distinction between the direct im- 
pact ionization and RESI emerges from our classical cal- 
culations [15] which allow both electrons to interact with 
the field and are, therefore, fully nonsequential. Evidence 
for the veracity of the proposed mechanism comes from 
two observations: The ionization delays are in harmony 



with experimental observations, and our treatment con- 
firms that the second ionization in RESI happens at the 
extrema of the laser field |5l llQ^16-18|. 

Classical ensemble methods [19|, |20| have been remark- 
ably successful in identifying direct impact ionization as 
well as RESI pathways that lead to double ionization 
and reproduce the experimental and computational ob- 
servations closely. Therefore we consider a two active 
electron atom with soft-Coulomb potentials subjected to 
an intense and short linearly polarized laser pulse in the 
dipole approximation 15L ll9M22| . The Hamiltonian is: 



H(xi,X2,Pi,P2,t) 
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+ (xi + X2) • e^Eofit) sincjt, 



(1) 



where x^ is the position vector of the zth electron in d di- 
mensions and Pi is its canonically conjugate momentum. 
The linearly polarized (along the x-direction with unit 
vector e^:) laser field is characterized by its amplitude £^0 
and has a wavelength of 780 nm or 460 nm {uj = 0.0584 
or 0.1 a.u. respectively) with a shape f {t) consisting 
of two-cycle linear ramp-up and six laser cycle constant 
plateau. The constants a and b are the electron-nucleus 
and electron-electron softening parameters respectively. 
For the commit at ions considered in this paper, we choose 
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15|, [19-22], even though qualitatively similar 



results are observed with different softening parameters. 

To begin with, we investigate the dynamics when both 
electrons are confined to a single dimension along the 
axis of polarization [23], i.e., d = 1. The motivation 
for studying the one-dimensional model is twofold: First, 
the laser field drives the dynamics along the polarization 
axis by which ionization is naturally expected to hap- 
pen. Second, we find that this simplified dynamics forms 
the skeleton of higher-dimensional dynamics, which, how- 
ever, differs from the single-dimensional dynamics in sig- 
nificant ways, to be specified below. It should be noted 
that, three-dimensional calculations give similar results 
to two-dimensional ones due to the cylindrical symmetry 
of the problem around the polarization axis. 

RESI events are described by two-electron trajectories 
for which, after a recollision, one electron remains bound 
to the nucleus for a long time before ionizing without any 
additional recollision, i.e., without the other electron ever 
returning to the nucleus. We use an energy criterion to 
designate an electron as ionized or not [15, 24]. A recol- 
lision is said to have occurred whenever the distance be- 
tween the two electrons is smaller than some threshold. 
We consider the NSDI trajectories which double-ionize 
after a last recollision. Among them, we discriminate RE- 
SIs by a time delay of at least two laser cycles between 
the last recollision and the ionization of the remaining 
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FIG. 2: Stroboscopic plots at the maxima of the laser field, 
for / = 3 X 10^^ W • cm~^ and 780 nm during the plateau 
of the laser. Upper panel: detected RESI trajectories, after 
the last recollision, for Hamiltonian ([1]). Initial conditions are 
chosen as in Refs. [13, [23. For each RESI, the first point on 
the section is plotted in light blue while the following ones 
are in dark blue. Lower panel: stable (light orange dots) 
and unstable (dark blue dots) manifolds of the periodic orbit 
O12 of Hamiltonian ([2|). The position of O12 is indicated 
with a black diamond (see the arrow). The pink area in the 
panels represents the part of phase space from which the inner 
electron does not ionize. 



ion. Qualitatively similar results are obtained with dif- 
ferent (shorter or longer) time delays. In Fig. [2] (upper 
panel), we display stroboscopic coordinate- momentum 
plots (taken at the maxima of the field) of detected RE- 
SIs. In what follows, we explain the rationale for these 
swirling patterns and connect them to full-dimensional 
calculations. 

RESIs are best understood using a single electron 
model: A recollision has excited the inner electron while 
the recolliding electron remains ionized [15| , which allows 
us to neglect the electron-electron interaction in Hamilto- 
nian ([1]). The resulting reduced-dimensional Hamiltonian 



reads 
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where (f)o denotes the phase at which the last recohision 
occurs. 

The dynamics given by Hamiltonian (J2j) turns out to 
be the key to the patterns formed by the RESI trajecto- 
ries. Previous studies on Hamiltonian (J2j) [15|, [2^ with 
one spatial dimension have identified two qualitatively 
different kinds of dynamics for the electron driven by the 
field: In the competition between Coulomb attraction to 
the nucleus and the laser excitation, either the latter pre- 
vails, and the electron is quickly ionized; or the Coulomb 
attraction manages to maintain the electron trapped near 
the core [25| . Two areas in phase space emerge from this 
distinction: A bound region, close to the nucleus (pink 
area in Fig. [2]), where the electron is trapped by the nu- 
cleus and cannot ionize; and an unbound region, further 
away (white area in the same panels), where the electron 
is quickly ionized by the laser. 

A more detailed study shows that the behavior in the 
unbound region is more complex than anticipated. This 
is readily apparent in the patterns seen in the unbound 
region in Fig. O A thin transition layer in the unbound 
region, and located in the area where the Coulomb at- 
traction and the laser excitation compete equally, is re- 
sponsible for RESI. In practice, we show that this transi- 
tion region is organized by the main resonances between 
the free field dynamics [Eq = in Hamiltonian (J2])] and 
the laser (the electron revolves exactly n times around 
the nucleus in one laser cycle, and we refer them as l:n 
resonances). These resonances give birth to periodic or- 
bits, among which at least one, denoted On, is unsta- 
ble (hyperbolic). Other periodic orbits in the vicinity of 
the bound region merely influence the fine details of the 
chaotic structure. 

Because a periodic orbit corresponds to a recurrent mo- 
tion, it does not lead to ionization. However, neighbor- 
ing trajectories may do so after some time, particularly 
if the periodic orbit is weakly unstable, as it is the case 
for selectes (9nS. The pathways by which the electrons 
approach or leave the core are quantified by the so-called 
stable and unstable manifolds [26] around an unstable 
periodic orbit. In Fig. [2] (lower panel) we draw the sta- 
ble and unstable manifolds of On (light orange and dark 
blue points respectively) forn = 12. Note the strong sim- 
ilarity between the unstable manifold and the Poincare 
sections of RESI trajectories for the two-electron Hamil- 
tonian ([T]) (upper panel). This similarity confirms the 
key role played by this one-electron unstable manifold in 
the RESI process. In addition, we see that the stable 
and unstable manifolds intersect an infinite number of 
times, a characteristic feature of a chaotic dynamics [27|. 
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FIG. 3: Projection of the unstable manifold of the periodic 
orbit Oi2 of Hamiltonian Q for / = 3 x 10^^ W • cm"^ and 
780 nm in the (x, 0) plane. The black curve is a projection of 
the periodic orbit in the plane (x,(/)), and the black squares 
indicate the position of the saddle point at the maximum of 
the field. 



The overlap between the stable and unstable manifolds 
of On forms a "sticky" region ^] that traps trajectories 
for some time before ionizing. To better understand how 
RESI electrons find their way to ionization, we display a 
projection of the unstable manifold of On in the position- 
phase of the laser plane {x , (j) = cot -\- (j)o) , together with 
a projection of On in Fig. [3l Two main branches depart 
from the central region near x = when the laser phase is 
(j) = 7r/2 and (j) = 37r/2, i.e., near the extrema of the elec- 
tric field, showing that RESI takes place approximately 
at the extrema of the electric field [5i, llTJ. 

The organizing role of resonant periodic orbits can also 
be found in the relative proportion of RESI to double ion- 
ization yields: In Fig. IH we compute the RESI yields as 
the intensity is varied. It shows oscillations which are cor- 
related with the intensity range where the l:n resonance 
regulates the delayed ionization. The bound region (pink 
area in Fig. [2]) shrinks with increasing laser intensity, and 
higher-order resonant periodic orbits are drawn to the un- 
bound region. It means that as On becomes too unstable, 
the next resonance (associated with periodic orbit On+i) 
is at play. Note the correlation between the oscillations 
and the stability index [26] of the resonant periodic orbits 
identified from the one-electron model (|2]). The stability 
index measures the typical time this orbit is expected to 
influence the neighboring dynamics: The larger the sta- 
bility index, the sooner a neighboring trajectory diverges 
from it. Recent experimental results [28], with elliptic 
polarization, have revealed oscillations in the parallel to 
anti-parallel double ionization yields analogous to the os- 
cillations observed in the RESI yields of Fig. |4j 

Adding a dimension to the one-dimensional calcula- 
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FIG. 4: Normalized RESI yields for 460 nm wavelength laser 
(markers, left hand ^/-axis). Red squares (diamonds) label 
one (two) dimensional one-electron simulations (normalized 
to the number of non- ionized trajectories). Green dots label 
two-dimensional two-electron simulations (normalized to the 
number of double ionization) . Curves label the linear stabil- 
ity index [2y] of the resonant periodic orbits On of Hamilto- 
nian (|2) (curves, right hand y-axis) as functions of the laser 
intensity /. Continuous (broken) curves refer to odd (even) 
l:n resonances. As intensity increases, the order of the reso- 
nance goes from n = 3 to n = 9 as indicated on the curves. 



tions provides more scope for the electron dynamics and 
makes the dynamics more complex. In particular, due 
to the extra dimension, the dynamics close to the nu- 
cleus region becomes partially chaotic both for one- and 
two-electron Hamiltonians (|2]) and ([1]) as it is already 
seen in Fig. [H It becomes harder to define a bound 
and unbound regions for the inner electron. More im- 
portantly, the aforementioned l:n resonances generally 
enhance (linear) instability in the transverse direction 
while not affecting the organization of the dynamics in 
the (polarization) symmetry subspace: At a given inten- 
sity more than one resonance drives the RESI dynamics 
in the sense that their unstable directions participate to 
the enhancement of delayed ionization. A direct conse- 
quence is the smoothing of the oscillations in Fig. 2] for 
460 nm wavelength. Increasing the wavelength to 780 nm 
gives birth to more resonances with a more dense tangle 
of unstable manifolds so that the oscillations are com- 
pletely washed out in two-dimensional models whether 
they are one- or two-electron models, as observed from 
our classical calculations. 

In Fig. H] we notice that the RESI yields are higher for 
a two-dimensional model than for one dimension. We at- 
tribute these higher yields to the aforementioned chaotic 
dynamics coming from an increase of instability associ- 
ated with additional resonances in the full-dimensional 
models. It underscores once more the pivotal role played 



by the inner electron dynamics [Hamiltonian (J2]), i.e., 
without recollision] in these delayed ionizations. 
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